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Liesegang  Rings  and  a  Theory  of  Fast 
Reaction  and  Slow  Diffusion 

Joseph  B.  Keller 


1.  INTRODUCTION. 


A  chemically  reacting  system  usually  proceeds  mono- 
tonically  in  time  toward  a  state  of  completed  reaction ,  a 
state  of  equilibrium,  or  a  steady  state.  The  ultimate  state 
is  generally  homogeneous  in  space,  as  a  consequence  of  diffu¬ 
sion,  even  though  the  initial  state  may  be  inhomogeneous. 
Therefore  it  is  surprising  that  certain  chemical  systems  do 
not  behave  in  this  way.  Instead  some  of  them  exhibit  temporal 
oscillations  which  do  not  decay  in  time,  but  which  become 
temporally  periodic.  Other  systems,  or  the  same  systems  with 
different  parameter  values,  perform  oscillations  which  persist 
but  remain  chaotic  rather  than  becoming  periodic.  Still 
others  reach  steady  states  which  are  not  homogeneous  in  space, 
but  which  are  spatially  periodic  or  quasi-periodic.  There  are 
also  systems  which  continue  to  oscillate  in  time  and  to  be 
quasi-periodic  in  space. 

Examples  of  such  unusual  behavior  of  chemically  reacting 
systems  have  been  known  for  a  long  time.  One  of  them  was  dis¬ 
covered  by  R.  E.  Liesegang  [  1)  in  1896  and  studied  by  him  for 
many  years  thereafter.  The  spatial  pattern  which  occurs  in 
this  system  consists  of  a  family  of  concentric  circular  rings 
which  are  called  Liesegang  rings  in  his  honor.  It  is  this 
system  and  similar  ones  which  we  shall  examine.  In  order  to 
analyze  them,  first  we  shall  have  to  describe  in  detail  the 
physical  and  chemical  processes  which  occur.  This  we  shall  do 
following  the  ideas  of  W.  Ostwald  [2]  who,  in  1897,  outlined 
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an  explanation  involving  precipitation  and  supersaturation. 

We  shall  extend  that  out  line  to  a  complete  description  of  the 
various  processes.  Then  we  shall  convert  this  verbal  descrip¬ 
tion  into  a  mathematical  theory,  and  formulate  an  appropriate 
mathematical  problem  corresponding  to  it.  All  of  these  con¬ 
siderations,  and  the  subsequent  analysis  of  the  problem,  are 
the  joint  work  of  S.  I.  Rubinow  and  the  author,  and  are  pre¬ 
sented  in  greater  detail  in  their  paper  [3]. 

In  order  to  treat  the  mathematical  problem  we  shall  make 
certain  appropriate  simplifications.  The  most  important  of 
these  is  a  consequence  of  the  fact  that  the  reaction  process 
is  fast  compared  to  the  diffusion  process.  This  has  led  us  to 
formulate  and  develop  a  general  theory  of  reaction  and  diffu¬ 
sion  processes  involving  fast  reaction  and  slow  diffusion.  We 
shall  describe  this  theory,  which  should  have  many  other 
applications,  and  use  it  to  simplify  the  present  problem.  It 
was  worked  out  by  P.  S.  Hagan  and  the  author,  and  is  contained 
in  their  paper  [4]  together  with  the  related  theory  of  some 
fast  and  some  slow  reactions. 

The  subsequent  analysis  of  the  simplified  problem  is 
carried  out  in  [3]  employing  approximations  like  those  in  the 
related  work  of  C.  Wagner  [5],  S.  Prager  [6]  and  Zeldovich, 
Barenblatt  and  Salganik  [7].  Some  predictions  of  the  theory 
are  compared  in  [3]  with  certain  relevant  experimental  results . 
More  detailed  comparison  between  theory  and  experiment  will  be 
possible  when  the  equations  of  the  theory  are  solved 
numerically.  J.-M.  Vanden-Broeck  and  the  author  [8]  are 
solving  them  by  using  the  method  of  finite  differences  to  re¬ 
place  the  partial  differential  equations  by  a  finite  set  of 
algebraic  equation's.  These  equations,  together  with  the 
initial,  boundary,  and  jump  conditions*  are  being  solved  by  an 
appropriate  iterative  method. 

Before  turning  to  the  study  of  Liesegang  rings,  we  shall 
comment  on  the  periodic  non-monotonic  behavior  of  chemical 
systems  in  general.  Let  us  first  consider  why  chemical  reac¬ 
tions  are  expected  to  proceed  monotonically  in  time.  The 
reason  is  that  the  rate  of  reaction  is  positive  as  long  as  the 
concentrations  of  the  reactants  are  positive  and  the  concen¬ 
tration  of  the  reaction  products  is  not  too  large.  Furthermore 
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there  is  no  mechanism#  such  as  inertia,  to  cause  the  reaction 
to  overshoot  the  equilibrium  state  or  the  state  of  completed 
reaction.  Therefore  it  is  to  be  expected  that  a  reaction  will 
proceed  in  one  direction  at  a  decreasing  rate.  However  when 
several  reactions  occur,  there  is  no  longer  any  justification 
for  this  belief.  The  presence  of  other  processes,  such  as 
flow  and  diffusion,  can  also  lead  to  destabilization  of  uni¬ 
form  steady  states  and  the  occurrence  of  temporal  or  spatial 
periodicity. 

The  interest  in  periodic  and  non-uniform  states  of  reac¬ 
ting  systems  has  been  generated  by  scientific  curiosity,  by 
engineering  necessity,  and  by  potential  applications  in 
biology.  The  engineering  reason  is  the  need  to  control  and 
eliminate  oscillations  which  arise  in  stirred  tank  and  other 
continuous  flow  reactors.  The  biological  reasons  include  the 
possible  explanation  of  biological  clocks  and  a  proposed  ex¬ 
planation  of  morphogenesis.  The  latter  of  these,  made  by 
A.  M.  Turing  in  1952,  is  that  chemical  reaction  and  diffusion 
are  responsible  for  the  creation  of  form  in  biological  systems 
through  the  development  of  spatially  nonuniform  steady  states. 
All  of  these  reasons  have  led  to  extensive  experimental  and 
theoretical  investigations  of  the  non-monotonic  and  non-uniform 
behavior  of  chemical  systems. 

2.  THE  MECHANISM  OF  LIESEGANG  RING  FORMATION » 

Liesegang,  while  experimenting  with  photographic 
materials,  prepared  a  gel  containing  potassium  dichromate. 

This  gel  was  in  the  form  of  a  thin  layer  on  a  glass  plate.  He 
placed  a  drop  of  silver  nitrate  solution  on  the  gel  and  found, 
after  some  time,  that  concentric  circular  rings  of  silver 
chromate  formed  in  the  gel.  The  process  whereby  the  silver 
chromate  formed  into  visible  solid  particles  is  precipitation, 
so  this  phenomenon  has  sometimes  been  called  “periodic" 
precipitation.  However  it  is  not  periodic  in  either  space  or 
time,  so  that  is  a  misnomer.  We  prefer  to  call  it  recurrent 
precipitation,  since  the  rings  form  one  after  another  with  in¬ 
creasing  spatial  and  temporal  intervals  between  successive 
rings . 
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To  explain  this  phenomenon  we  suppose  that  silver  ions 
from  the  drop  go  into  solution  in  the  gel  and  diffuse 
radially  outward.  Within  the  gel  they  encounter  chromate  ions 
and  react  with  them  to  form  silver  dichromate.  As  the 
chromate  ions  near  the  drop  are  consumed  by  the  reaction,  more 
chromate  ions  diffuse  radially  inward.  Thus  the  two  reactants 
diffuse  toward  one  another  and,  upon  coming  together, 
continually  produce  silver  dichromate.  In  this  way  an 
expanding  circular  region  of  silver  dichromate  is  produced. 
This  is  the  monotonic  behavior  to  be  expected  in  chemically 
reacting  systems.  The  question  then  arises  of  why  only 
separated  rings  of  silver  dichromate  are  visible  at  the  end  of 
the  experiment? 

To  answer  this  question  we  assume  that  the  silver 
dichromate  is  not  visible  until  it  precipitates  out  of  solu¬ 
tion  and  forms  solid  particles.  Now  it  is  well  known  that  a 
dissolved  substance  cannot  precipitate  until  its  concentration 
c  exceeds  the  saturation  concentration  cs.  Then  the  rate  of 
precipitation  is  proportional  to  c  -  cs,  provided  that  some  of 
the  solid  is  present  upon  which  it  can  precipitate.  However 
if  none  of  the  solid  is  present  the  concentration  c  may  have 
to  reach  a  higher  value  c*  >  cs  before  precipitation  can  start. 
Therefore,  following  Ostwald,  we  assume  that  precipitation 
starts  only  when  c  >_  c*.  After  it  has  started,  we  assume  that 
it  continues  as  long  as  c  ^  cs. 

In  view  of  these  assumptions,  precipitation  will  begin  in 
the  gel  at  the  location  of  the  drop  of  silver  nitrate  solution , 
provided  that  the  reactants  are  present  in  sufficiently  high 
concentrations  to  produce  silver  dichromate  with  a  concentra¬ 
tion  c  ^  c*.  Precipitation  will  continue  within  a  circle  of 
increasing  radius  surrounding  the  drop  until  c  falls  below  c* 
at  the  outer  boundary  of  this  circle.  Then  the  precipitation 
zone  will  stop  growing. 

The  reason  why  c  may  fall  below  c*  is  that  the  chromate 
ions  in  the  gel  are  depleted  by  reaction  and  by  diffusion 
toward  the  silver  nitrate  drop,  while  the  concentration  of 
silver  ions  is  low  far  from  the  drop.  As  time  goes  on  this 
latter  concentration  increases  everywhere.  A  time  may  be 
reached  when  this  silver  ion  concentration  and  that  of  the 
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chromate  ions  is  high  enough  somewhere  in  the  gel  to  produce 
silver  dichromate  with  a  concentration  c  ^  c*.  Then  a  second 
zone  of  precipitation  will  start  to  form  at  that  place.  After 
a  while  this  process  will  be  repeated,  forming  a  third  zone, 
and  so  on.  The  resulting  zones  are  the  Liesegang  rings. 

Since  Liesegang' s  first  experiment,  the  same  phenomenon 
has  been  found  to  occur  in  many  other  reactions.  Furthermore, 
it  has  also  been  shown  to  occur  in  test  tubes  and  capillary 
tubes,  with  the  formation  of  planar  bands  at  various  positions 
along  the  tube  ,^?tead  of  rings.  It  is  this  one  dimensional 
case  which  we  shalx  treat. 

3.  A  MATHEMATICAL  FORMULATION  OF  THE  THEORY. 

To  formulate  the  preceding  theory  mathematically,  we 
denote  by  a,  b,  c  and  d  the  molar  concentrations  of  the  four 
substances  A,  B,  C  and  D.  Here  A  is  analogous  to  the  silver 
ions,  B  to  the  chromate  ions,  C  to  the  dissolved  silver 
dichromate  and  D  to  the  solid  or  precipitated  silver  dichro¬ 
mate.  We  assume  that  moles  of  A  can  combine  with  vfi  moles 
of  B  to  form  vc  moles  of  C,  and  that  this  reaction  is  rever¬ 
sible.  We  also  assume  that  C  can  precipitate  to  form  D,  but 
that  D  cannot  dissolve  to  form  C.  In  addition  we  suppose  that 
A,  B  and  C  can  diffuse  with  diffusion  coefficients  D^,  Dfi  and 
Dc,  but  that  D  cannot  diffuse.  Then  the  four  concentrations, 
which  depend  upon  the  time  t  and  the  coordinate  x  along  the 
tube,  satisfy  the  following  reaction-diffusion-precipitation 
equations  within  the  tube  x  >  0: 

at  "  Vxx  •  Y  '  i3-1] 

bt  =  °BbXX  -  VBr  '  (3-2) 

Ct  =  DC°xx  +  vcr  -  P(c'd)  '  (3-3) 

dt  *  p (c,d)  .  (3.4) 


Here  r  is  the  reaction  rate  and  p  is  the  precipitation  rate. 

We  shall  assume  that  r  is  given  by  the  law  of  mass  action 
with  rate  constant  k+  for  the  forward  reaction  and  k_  for  the 
backward  reaction.  Thus 
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r (a,b,c) 


v  v 

k+a  Ab  B 


k  c 


(3.5) 


The  precipitation  rate  was  described  in  Section  2,  and  that 
description  leads  to  the  following  expression  for  p: 


p  (c  ,d)  =  0  ,  if  c  <  c*  and  d  *  0  , 

=  q(c-cs)+  ,  if  c  >  c*  or  d  >  0  . 


(3.6) 


Here  q  is  a  rate  constant  and  the  subscript  +  on  c  -  cs  de¬ 
notes  the  positive  part  of  c  *  cs. 

Since  we  have  assumed  that  the  reactions  occur  in  the 
tube  occupying  the  region  x  >  0,  we  must  specify  the  initial 
conditions  within  the  tube  and  boundary  conditions  at  the  end 
x  -  0.  To  correspond  with  Liesegang's  experiment,  we  assume 
that  only  B  is  present  initially.  Thus  we  require  that 


a  -  c  =  d  =  0  ,  b  s  bQ  at  t  =  0  ,  x  >  0  .  (3.7) 

At  the  endpoint  x  =  0  we  assume  that  A  is  kept  at  the  concen¬ 
tration  aQ  while  B  and  C  cannot  leave  the  tube  there.  Thus 

a  =  aQ  ,  bx=cx=0atx=0,  t>0.  (3.8) 

In  addition  to  these  conditions,  we  require  that  a,  b,  c  and 
their  x  derivatives  be  continuous  throughout  the  tube  for 
t  >  0: 

a,  b,  c,  ax,  bx,  cx  continuous  for  x  >  0,  t  >  0  .  (3.9) 

The  preceding  equations  constitute  our  mathematical  for¬ 
mulation  of  the  theory  of  Liesegang  band  formation.  The  bands 
are  the  regions  within  which  d(x,t)  is  positive  after  a 
sufficiently  long  time.  Thus  the  bands  are  defined  by 

d(x,»)  >  0  .  (3.10) 


This  formulation  is  complete  in  the  sense  that  the  equa¬ 
tions  seem  to  have  a  unique  solution.  This  solution  deter¬ 
mines  exactly  where  the  bands  will  be,  the  thickness  of  each 
band,  the  time  when  each  band  starts  to  form,  the  concentra¬ 
tions  between  the  bands,  etc.  The  solution  will  also  show 
whether  or  not  bands  do  form,  and  it  will  yield  the  sets  of 


WfliWJWWJflW* 


LIESEGANG  RINGS;  THEORY  OF  FAST  REACTION  AND  SLOW  DIFFUSION  217 


parameter  values  for  which  they  form.  The  solution  can  be 
calculated  numerically,  but  we  have  not  calculated  it. 
Instead  we  have  first  simplified  the  equations  and  then 
analyzed  the  resulting  problem. 


4.  REDUCTION  TO  A  SIMPLER  PROBLEM. 

In  order  to  simplify  the  theory  presented  in  Section  3, 
we  shall  proceed  heuristically ,  making  assumptions  based  upon 
intuitive  chemical  considerations.  Later  on  we  shall  indicate 
how  the  same  simplified  theory  results  from  a  systematic 
application  of  the  theory  of  fast  reactions  and  slow  diffusion. 
Our  first  assumption  is  that  the  reaction  term  v^r  in  (3.1)  is 
negligible  compared  to  the  other  terms  in  that  equation.  This 
assumption  is  valid  when  the  concentration  of  B  is  sufficiently 
small.  It  yields 

a  =  Da  .  (4.1) 

t  A  xx 

Next  we  add  v_  times  (3.3)  to  v  times  0.1)  to  eliminate 

B  L 

r  and  obtain 


'•cbt  +  VBCt  "  VCDBbxx  +  VBDC°xx  -  VBP(C'd)  • 


(4.2) 


Then  we  assume  that  the  reaction  proceeds  so  rapidly  that  it 
is  in  equilibrium,  which  yields 

r (a, b,c)  =  0  .  (4.3) 

We  shall  use  (4.2)  and  (4.3)  instead  of  (3.2)  and  (3.3).  But 
then  we  must  replace  the  two  boundary  conditions  =  cx  “  0 
by  the  single  condition  that  B,  both  free  and  contained  in 
C,  cannot  escape  through  the  boundary: 

VCDBbx  +  VB°CCx  "  0  at  x  =  0  ,  t  >  0  .  (4.4) 


With  these  modifications  the  problem  can  be  dealt  with  by 
first  solving  (4.1)  for  a(x,t)  with  a(x,0)  -  0  and  a(0,t)=ao. 
The  solution  is 

a(x,t)  =  aQ  erfc[x/(4DAt) ]  .  (4.5) 


Then  (4.3)  can  be  solved  for  b  in  terms  of  a  and  c  and  the  re¬ 
sult  can  be  used  to  eliminate  b  from  (4.2)  and  (4.4).  In  this 
way  (4.2)  becomes  an  equation  for  c  in  which  d  seems  to  occur. 
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However  both  in  the  region  where  p  =  0  and  in  the  region  where 
p  >  0,  d  does  not  occur  in  the  equation  for  c.  Thus  the 
entire  problem  can  be  reduced  to  one  of  determining  c  in  each 
of  these  regions.  The  boundary  between  the  two  regions  con¬ 
sists  of  the  curve  on  which  c  =  c*  together  with  the  upper 
halves  of  the  tangents  to  this  curve  where  it  becomes 
vertical.  Once  c  is  found,  d  can  be  calculated  from  (3.4). 

The  reduction  just  described  is  carried  out  in  detail  in 
the  work  of  Keller  and  Rubinow  f3],  where  the  resulting  equa¬ 
tions  are  analyzed  and  solved  approximately.  Therefore  we 
shall  not  describe  that  analysis  and  its  results.  Instead  we 
shall  present  the  theory  of  fast  reaction  and  slow  diffusion, 
and  show  how  it  '"an  be  used  to  derive  the  reduced  problem  in  a 
systematic  way. 

5.  FAST  REACTION  AND  SLOW  DIFFUSION. 

In  many  chemical  systems,  such  as  those  considered  in  the 
preceding  sections,  both  chemical  reactions  and  diffusion 
occur,  but  at  quite  different  rates.  Often  the  reactions  pro¬ 
ceed  much  more  rapidly  than  the  diffusion.  In  such  cases  it 
is  possible  to  simplify  the  analysis  of  the  process  by  taking 
advantage  of  this  difference  in  rates.  One  way  of  doing  so  is 
to  assume  that  the  reactions  have  reached  equilibrium  or 
completion  at  each  point  of  space.  Then  the  state  at  each 
point  is  constrained  to  satisfy  the  condition  that  the  reac¬ 
tion  rates  vanish  there.  When  this  constraint  is  adjoined  to 
the  reaction-diffusion  equations  governing  the  system, an  over¬ 
determined  set  of  equations  results.  Chemical  and  physical 
reasoning  can  then  be  used  to  delete  some  of  the  equations  and 
thereby  make  the  system  determined. 

Instead  of  using  this  intuitive  procedure,  it  is  possible 
to  proceed  in  a  systematic  mathematical  manner  which  we  shall 
now  describe.  The  basic  idea  is  to  introduce  into  the  equa¬ 
tions  a  small  parameter  e  which  is  the  ratio  of  the  time  scale 
for  reaction  to  the  time  scale  for  diffusion.  Then  the 
asymptotic  expansion  of  the  solution  is  sought  for  e  near  zero. 
This  expansion  can  be  constructed  either  by  the  "two  time" 
method  or  by  the  method  of  matched  asymptotic  expansions. 
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Both  of  these  methods  are  utilized  in  [4],  Here  we  shall 
indicate  how  to  obtain  the  desired  results  via  the  latter 
method. 

Let  us  begin  by  considering  the  following  system  of 
reaction-diffusion  equations: 

ut(x,t,e)  =  f(u,x,t,et,e)  +  eDAu  ,  (5.1) 

u (x, 0 , e  )  =  g  (x)  .  (5.2) 

Here  u(x,t,e)  is  a  vector  of  concentrations#  temperature  and 
possibly  other  dependent  variables,  f (u,x, t # et,e )  is  a  vector 
of  reaction  rates,  heat  production  rate,  etc.  and  D  is  the 
diffusion  coefficient  matrix.  The  essential  feature  of  (5.1) 
is  that  the  diffusion  term  contains  the  factor  e,  which  is 
small.  We  have  also  permitted  f  to  depend  upon  all  the 
variables,  and  even  to  have  a  slow  time  dependence  via  the 
argument  et.  In  the  initial  condition  (5.2)  the  function  g(x) 
is  given. 

To  obtain  the  asymptotic  expansion  of  u  for  e  small,  we 
assume  first  that  it  is  of  the  form 

u(x,t,e)  -  u°(x,t)  +  eu^  (x,t)  +  •••  .  (5.3) 

We  shall  call  this  the  initial  layer  expansion  of  u.  We  now 
substitute  (5.3)  into  (5.1)  and  (5.2)  and  set  e  =  0  to  obtain 

u°(x,t)  =  f (u°,x,t,0,0)  ,  (5.4) 

u°(x,0)  =  g (x)  .  (5.5) 

These  equations  describe  the  evolution  in  time  of  u°  at  the 
point  x.  Further  terms  in  (5.3)  can  be  found  by  considering 
the  coefficients  of  higher  powers  of  e  in  (5.1)  and  (5.2),  but 
we  shall  not  examine  them.  Instead  we  shall  consider  the 
behavior  of  u°(x,t)  as  t  +  ®. 

To  this  end,  we  assume  that  f (u^,x, t, 0, 0)  has  a  limit  as 
t  ■+■  00 ,  and  then  we  may  expect  u^(x,t)  to  have  a  limit  also. 

Let  us  denote  it  by  u^(x,°»).  Thus  we  assume  that 

lim  u°(x,t)  =  u0(x,°°)  . 

t>® 


(5.6) 
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Then  u° (x,°°)  must  satisfy  the  following  equilibrium  condition, 
which  results  from  (5.4): 

f [u°(x.®) fx,», 0,0]  =  0  .  (5.7) 

In  general  u0(x,°°)  is  not  uniquely  determined  by  (5.7),  but  it 
is  determined  by  the  initial  condition  (5.5)  and  the  reaction 
equation  (5.4).  It  is  this  limit  u°(x,°°)  of  u°(x,t)  which 
will  be  needed  to  determine  the  outer  expansion  of  u,  which  we 
consider  next. 

The  initial  layer  expansion  (5.3)  is  presumably  valid  for 
fixed  values  of  t  as  e  +  0.  In  order  to  find  the  behavior  of 
u  for  long  times  of  order  e  ^ ,  we  must  construct  another 
expansion,  called  the  outer  expansion  of  u.  In  order  to  con¬ 
struct  it  we  first  introduce  the  new  time  variable  t  =  et  and 
consider  u  to  be  a  function  of  t: 

u(x,t,e)  =  v(x,x,e)  ,  x  =  et  .  (5.8) 

Then  (5.1)  becomes  the  following  equation  for  v: 

evT(x,x,e)  =  f (v,x,e-1xfT , e )  +  eDAv  .  (5.9) 

Next  we  assume  that  v  has  the  expansion 

v(x,T,e)  =  v°(x,t)  +  ev^1^  (x,t)  +  .  (5.10) 

This  is  the  outer  expansion  of  u. 

Both  the  inner  expansion  (5.5)  and  the  outer  expansion 
(5.10)  represent  u.  If  there  is  a  common  region  of  validity 
of  the  two  expansions,  as  we  assume,  then  in  it  they  must  be 
asymptotic  to  one  another.  Thus  we  have  in  this  region 

u°(x,t)  +  eu^(x,t)  +  •••  **  v°(x,et)  +  ev^(x,et) 

+  •••  .  (5.11) 

-1/2 

Suppose  that  (5.11)  holds  for  t  =  e  '  .  Then  by  choosing 

this  value  for  t  and  letting  e  ->  0  in  (5.11)  we  obtain 

u°(x,»)  =  v° (x, 0)  .  (5.12) 

This  procedure  for  deriving  (5.12)  is  called  "matching"  the 
initial  layer  expansion  and  the  outer  expansion.  The  result 
(5.12)  yields  the  initial  value  of  v^  terms  of  the  "final" 
value  of  u  . 
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We  next  substitute  (5.1)  into  (5.9)  and  set  e  =  0  to  get 
f [v°(x,t) 0]  =  0  .  (5.13) 

This  equation  shows  that  must  lie  on  the  equilibrium  sur¬ 
face  (or  manifold)  defined  by  f  =  0,  but  it  does  not  determine 
where  is  on  this  surface  unless  the  surface  is  a  single 
point.  From  (5.12)  we  know  that  at  t  =  0,  v°  starts  at  the 
point  u^(x,°°).  To  find  how  it  evolves  as  t  increases,  we 
consider  the  terms  of  order  ^  in  (5.9)  ,  which  yield  the 
following  linear  equation  .or  v^^: 

fy[  v°  (x,  t  )  ,  x,°° ,  r  ,  0  1  v  ^  ^  =  v°  -  DAv°  -  f£  .  (5.14) 

The  gradient  matrix  f  in  (5.14),  which  is  the  coefficient 
(1)  v 

of  v  ,  will  generally  be  a  singular  matrix.  Let  us  suppose 
that  it  has  k  left  null-vectors  .  Then  (5.14)  is 

(1)  IK 

solvable  for  v  only  if  the  right  side  satisfies  k 

solvability  conditions.  They  can  be  obtained  by  scalar  multi¬ 
plication  of  (5.14)  on  the  left  by  each  left  null-vector. 

This  yields  the  k  conditions 

£.*{v^(x,t)  -  DAv^(x,t)  -  f  [v° (x,x ) ,x,~,t , 0) }  =  0  , 

J  T  G 

j  =  l,...,k  .  (5.15) 

This  is  a  system  of  k  equations  that  determine  how  v°(x,x) 
diffuses  on  the  equilibrium  surface  upon  which  it  is  con¬ 
strained  to  lie  by  (5.13).  The  initial  value  of  v^1  is  given 
by  (5.12).  Thus  (5.12),  (5.13)  and  (5.15)  determine  the  slow 

diffusion  of  v^(x,t),  the  leading  term  in  the  outer  expansion 
of  u(x,e  ^T,e).  The  fast  reaction  is  governed  by  (5.4)  and 
(5.5),  which  determine  u^(x,t),  the  leading  term  in  the  initial 
layer  expansion  of  u(x,t,e) . 

The  system  of  equations  (5.13)  and  (5.15),  together  with 
the  initial  condition  (5.12),  are  the  results  of  our  analysis. 
We  shall  not  examine  them  further  here,  nor  shall  we  obtain 
more  terms  in  the  outer  expansion.  These  and  other  matters  are 
considered  in  [4], 
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6.  APPLICATION  TO  THE  THEORY  OF  LIESEGANG  BANDS. 

We  shall  now  apply  the  theory  of  Section  5  to  the  equa¬ 
tions  of  Section  3,  which  constitute  our  formulation  of  the 
theory  of  Liesegang  bands.  To  do  so  we  first  rewrite  these 
equations,  introducing  explicitly  the  small  parameter  e  to  in¬ 
dicate  the  relative  size  of  each  term.  In  (3.1)  for  example, 
the  diffusion  term  is  small  of  order  £,  say,  but  the  reaction 

term  is  supposed  to  be  smaller.  Therefore  we  shall  assume  it 
2 

to  be  of  order  e  .  Thus  we  write  (3.1)  in  the  form 


2 

=  eD,a  -  e  v  r 
t  A  xx  A 


(6.1) 


Similarly  we  write  (3.2)  and  (3.3)  in  the  form 


b.  =  eD  b  -  v„r 
t  B  xx  B 


(6.2) 


=  eD^c  +  v_r  -  ep  . 
t  C  xx  C 


(6.3) 


The  system  (6.1) -(6. 3)  is  of  the  form  (5.1)  with  the 
following  identifications: 


2  ) 

a 

“e  VAr 

b 

,  f  = 

-VBr 

c 

,v„r  -  ep. 

D  - 


db  0 


(6.4) 


Thus  the  leading  term  of  the  initial  layer  expansion, 
u°  =  (a°,b°,c°) ,  satisfies  (5.4)  which  becomes 


0  , 

-vBr(a°,b°,c0)  , 


,0,00, 

ct  “  Vcr  a  'b  'c  >  • 


(6.5) 

(6.6) 

(6.7) 


From  (3.7)  the  corresponding  initial  conditions  are 


a°(x,0)  =  c°(x,0)  =  0  ,  b°(x,0) 


bQ  ,  x  >  0 


(6.8) 


■  V-.  .-*v.  ■  *.- 
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With  these  initial  conditions  and  r  given  by  (3.  S'),  the  solu¬ 
tion  of  (6-5)— (6.7)  is 

a°(x,t)  =  c°(x,t)  =  0  ,  b°(x,t)  =  bp  ,  x  >  0  .  (6.9) 

Thus  nothing  changes  in  the  initial  layer,  to  this  order. 

To  treat  the  outer  expansion  we  shall  write 


u(x,t,e)  =  v (x,t , e ) 


a  ( x ,  x  ,  e.) 
6 (x,t ,g) 
Y  ( x ,  t  ,  g  ) 


(6.10) 


Then  by  using  (6.4)  for  f  in  (5.13)  we  get  the  single  equation 


r[a°(x,x) ,6°(x,t) ,y°(x,t) ]  =  0  .  (6.11) 

Thus  the  leading  term  in  the  outer  expansion  satisfies  (4.3), 
as  we  assumed  in  Section  4 . 

Next  from  (6.4)  we  compute  f  at  e  =  0,  and  obtain 


f  =  (f  ,fw,f  ) 
v  a  b  c 


0 

0 

0 

VBra 

‘VBrb 

“VBrc 

VCra 

vcrb 

VCrc 

(6.12) 


This  matrix  has  two  linearly  independent  left  null-vectors 

i1  =  (1,0,0)  ,  z2  =  (0,VC,VB)  .  (6.13) 


We  also  compute  f£  at  e  =  0,  which  is  given  by 
’  01 

fe  =  0  • 

.-P. 


(6.14) 


We  can  now  use  each  of  the  vectors  and  given  by 
(6.13)  in  (5.15),  together  with  f^  given  by  (6.14)  and  D  given 
by  (6.4).  By  using  l ^  in  (5.15),  we  obtain 


a0  =  Daa° 
t  A  xx 


(6.15) 


This  is  just  (4.1)  of  Section  4,  with  a  replaced  by  a  and  t 
replaced  by  t*  Then  by  using  l2  (5.15)  we  get 

0  _,_.0 


VceT  +  VBY?  =  VcVxx  +  VBDCYxx  - 


(6.16) 
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This  is  exactly  (4.2)  with  b,  c  and  t  replaced  by  B^,  y°  and  i 
respectively. 

We  have  now  shown  how  the  simplified  equations  of 
Section  4  follow  from  the  equations  for  the  outer  asymptotic 
expansion  derived  in  Section  5.  The  corresponding  initial 
conditions  are  obtained  by  using  (6.9)  in  (5.12),  and  they  are 
just  those  used  in  Section  4.  An  extension  of  the  theory  of 
Section  5  is  needed  to  derive  the  boundary  conditions  (4.4) 
and  a0(0,i)  =  a  ,  but  we  shall  not  present  it  here. 

This  completes  the  application  of  the  theory  of  fast 
reaction  and  slow  diffusion  to  the  equations  of  Section  3, 
when  they  are  written  in  the  form  (6.1)  —  (6*3) .  Of  course  that 
way  of  writing  the  equations  depends  upon  the  parameters  in 
the  problem.  They  must  be  such  that  the  terms  in  these 
equations  have  the  indicated  relative  magnitudes. 
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